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Abstract 

We reconsider constraints on the eigenvalue density of the Dirac operator in the chiral symmetric 
phase of 2 flavor QCD at finite temperature. To avoid possible ultra-violet (UV) divergences, we 
work on a lattice, employing the overlap Dirac operator, which ensures the exact "chiral" symmetry 
at finite lattice spacings. Studying multi-point correlation functions in various channels and taking 



their thermodynamical limit (and then taking the chiral limit), we obtain stronger constraints than 
those found in the previous studies: not only the eigenvalue density at the origin but also its first 
and second derivatives vanish in the chiral limit of 2 flavor QCD. In addition we show that the axial 
U(l) anomaly becomes invisible in susceptibilities of scalar and pseudo scalar mesons, suggesting 
that the 2nd order chiral phase transition with the 0(4) scaling is not realized in 2 flavor QCD. 
Possible lattice artifacts when non-chiral lattice Dirac operator is employed are briefly discussed. 
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I. INTRODUCTION 



The QCD Lagrangian with Nf massless quarks is invariant under SU(Nf)i, x SU(Nf)n x 
U(l)y x U(1)a chiral rotations. This symmetry, however, is broken in two different ways: 
SU(Nf) L x SU(Nf) R part is spontaneously broken to SU(Nf) v in the QCD vacuum, while 
U(1)a is broken explicitly at quantum level by the anomaly. 

At a finite temperature T, it is widely believed that the SU(Nf)i, x SU(Nf)n chiral sym- 
metry is recovered above the (critical) temperature T c ~ 150 MeV and plenty of evidences 
have been reported in the first principle calculations of lattice QCD. For the U(1)a part, 
however, it remains an open question if, how (much), and when the symmetry is restored. 
We only know that the U(1)a symmetry should be recovered in the T — > oo limit, where 
fermions eventually decouple as the lowest Matsubara frequency goes to infinity, so that the 
anomaly term cannot survive. 

In particular, the question whether the U(1)a symmetry is restored or not near T c is of 
phenomenological importance. For simplicity, let us consider the Nf = 2 case. As Pisarski 
and Wilczek [l| have discussed, the order of the phase transition may depend on the fate of 
the U(1)a symmetry: if it remains to be broken at T c , the chiral phase transition can be the 
second order, while it is likely to be the first order when the U(1)a symmetry is also restored. 
Furthermore the particle spectrum with presence or absence of the U(1)a symmetry is quite 
different j^j] . A connection between the restoration of U(1)a symmetry and the gap in the 
eigenvalue density of Dirac operator near the origin is also suggested [3|]. 

In principle, the fate of the U(1)a symmetry and related issues can be investigated by 
numerical lattice QCD simulations |4|, |5[. Such studies are, however, still not easy, since 
both chiral and thermodynamical (the infinite volume) limits are required. Currently, four 
simulations with different quark actions are on-going, but they have reported different re- 



sults. Two of them 



, [7| have reported that the eigenvalue density of the Dirac operator has 
no gap at the origin and its quark mass scaling is consistent with the U(1)a broken scenario. 
Another group [8[ has also reported no gap at the origin but they have found that small 
eigenmodes, which mainly contribute to U(1)a breaking correlation functions, are localized 
and uncorrelated, suggesting that their contribution to the correlation functions is negligi- 



ble. A simulation with overlap quarks 



9[, however, has reported an existence of a gap in 



the Dirac eigenvalue density and a degeneracy of pion and eta(-prime) meson correlators, 
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which suggests the recovery of U(1)a symmetry. 

In this work, we address these problems again on a lattice, but using an analytic method. 
For simplicity, we concentrate on the Nf = 2 case in this paper. We employ the overlap 



Dirac operator 
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11], which ensures the exact SU(2)l x «S77(2)_r symmetry 12j through 



measure 



the Ginsparg- Wilson relation [13] but the U(1)a is (correctly) broken by the fermonic 



14|. By using the spectral decomposition of the multi-point correlation functions, 
and assuming the restoration of the SU(2)l x SU(2)r symmetry, we investigate if there are 
new constraints on the Dirac eigenvalue density in addition to the manifest one implied by 



the well-known Banks-Casher relation 15]. We also investigate whether the effect of the 
U(1)a symmetry breaking disappear or not above T c . 

Since similar analytical investigations have been made in the previous studies, let us here 
revisit them and make what's new in this work clear. The first analysis based only on QCD 



was done by Cohen [16]. Assuming an absence of the zero-mode's contribution, he concluded 
that all the disconnected contributions of the two-point functions in the SU(2)l X SU(2)r 
symmetric phase disappear in the chiral limit. This means that the pion, sigma, delta, 
and eta(-prime) meson correlators are all identical, realizing the U(1)a symmetry. In this 
work, we include the zero-mode contribution explicitly to check whether his conclusion still 
survives or not. 



In fact, Lee and Hatsuda |l7j (see also a related work by Evans et al. 18]) claimed that 
the zero-mode's contribution does not vanish but keeps the disconnected contribution of the 
scalar channel non-zero: 



m— ¥0 



Mm {(q(x)q(x) q(y)q(y)) - (q(x)T q(x) q(y)T q(y)) 



[ d[A} u=±1 e- s ^det'[f+m] 2 x A^(x)^(x)^(y)^(y) 
= lim/ + 0(m),(l) 

m— >0 Zj 

where q denotes the quark field, T 3 is the 3rd generator of SU(2). On the right-hand side 
(RHS), Z is the partition function of QCD, (i[A] iy=±1 means the gauge field integrals with a 
fixed topological charge v = ±1, Sym denotes the gauge part of the action, det'[^+ m] = 
det[p+ m]/m is the (continuum) fermion determinant with the quark mass m from which 
the zero-mode contribution is subtracted, and (J)q is an eigenfunction for the zero-mode at 
a given configuration A. 

The thermodynamical limit of Eq. ([1]) is, however, non-trivial and subtle, as was pointed 
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out by Cohen [3j. In fact, we find that the RHS of Eq. ([I]) is at least an 0(1 /V) quantity. 
Integrating Eq. (JT|) over x and then taking an average over y, (which should be grater than 
the original LHS), one immediately obtains 



since both = J d 4 x(j)Q (x)4>q (x) and C = lim m ^ (^i + Z_i)/m 2 Z (where Z±\ denotes the 
partition function in the topological sector of v = ±1) are finite. In fact, our work will show 
that not only Eq. (pQ) but also any contributions from zero modes of the Dirac operator are 
in general 0(1 /V) quantities, and thus disappear in the large volume limit. It is not difficult 
to intuitively understand our result. In the large volume limit that V — > oo , the number of 
the fermion modes contributing to the denominator Z increases (it is natural to assume it 
to be proportional to V ) while that of the numerator, where the bulk 0(V) contributions 
are canceled, is fixed to be 0(1). 

n 

Two years later, Cohen [3J discussed a constraint on the eigenvalue density of the QCD 
Dirac operator in the chiral limit. Relating the scalar one-point function and pseudoscalar 
two-point functions in the chiral symmetric phase, he concluded that the eigenvalue density 
near zero that 



must have a > 1. In this paper, we examine up to 4 point correlation functions in more 
various channels, and obtain a stronger constraint: a > 2. In the case of integer a, we 
believe that our constraint that a is equal to or larger than 3 should be the strongest, since 
we know a theory which has both a = 3 and unbroken SU(2)l x SU(2)r (and also U(1)a) 
chiral symmetries: 2-flavor massless free quarks. 

Although we perform no numerical analysis in our study, we would like to discuss possible 
artifacts in lattice QCD simulations. In our analysis, the fully recovered SU(2)l x SU(2)r 
symmetry is crucial. We discuss possible modifications to our conclusions due to discretiza- 
tion effects if a non-chiral quark action is employed in numerical simulations. We also 
comment on finite volume effects. 

Our paper is organized as follows. In Section HH we explain our set-up, what we observe, 
and what we assume. The constraints on the eigenvalue density with integer power at the 
origin are given in Section IIII1 In Section IIVt we then address a question on the fate of the 
U(1)a symmetry. In Section [V] we discuss possible systematic effects which may arise in 




V — y oo 



(2) 



p(X) ~ |A 



(3) 
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lattice QCD simulations. Section I VI I is devoted to a case where the eigenvalue density has 
a fractional power at the origin. Conclusion and discussion are given in Section IVHI Some 
useful formula and detailed calculations are collected in two appendices. 



II. LATTICE SETUP 

A. Spectral decomposition of the overlap fermion 

We consider iV/-flavor lattice QCD in a finite volume V, with the (anti-)periodic boundary 
condition in space(time). 

The quark part of the action is given by 

S F = a 4 \$D{A)iP + mi)F(D(A))i)\ (x), F(D) = 1 - ^D, (4) 

X ^ 

where ip = (ipx,^, ■ ■ ■ipN f ) T denotes the set of Nf fermion fields with the degenerate mass 



m, a is the lattice spacing, and D(A) is the overlap Dirac operator 10|, for a given gauge 
field A, 

D(A) = -L(l + - DwW-VB* V (5) 

Ra \ J(D W (A) - l/Ra)1(D w (A) - 1/Ra)) 



Here Dw{A) denotes the Wilson Dirac operator for the same gauge configuration A, and R 
is an arbitrary constant. We have omitted the identity matrix l^ fX ^ f for the flavor indices 
for simplicity. 

It is well-known that the overlap Dirac operator satisfies the 75 hermiticity, D(A)^ = 



7 5 D(A)7 5 , and the Ginsparg- Wilson (GW) relation [l3| 



^(^4)75 + 75^(^4) = aD(A)R 75 D(A). (6) 



12[ under 



With this relation, the action Eq. (j4]) in the m — > limit is exactly symmetric 
the lattice chiral rotation, 

8 a <iP(x) = z0T a 7 5 [(l - RaD(A))ifj}(x) J (7) 
6 a ij(x) = z^(x)T a75 , (8) 

where 9 is an infinitesimal real parameter, and T a denotes the generator of SU(Nf) for 
a = 1,2, ■ ■ ■ Nj — 1, and T (= ljV/xiV/) denotes that for U(1)a- 



We now consider eigenvalues and eigenfunctions of D(A): D(A)cj)^ = A^0^. The GW 
relation implies that 

X* + X* = aRXiXt (9) 

where A^, and its complex conjugate A^ are in general complex numbers and therefore 
(</)£)' D (Ay = ((frnV^n- Moreover, from the GW relation (jSJ) and its consequence © we 
have 

D(A)l^i = - ^^ 75^ = *S#f. (10) 

Since (A^ — A^) (0^,750^) = 0, eigenfunctions with complex eigenvalues can be ortho- 
normalized as (0^,0^) = (750^750™) = $nm, and (0^,7 5 0^) = 0. Here an inner product 
is defined as (/, g) = a A ^2 x p(x)g(x). For the real eigen values A^ = and A^ = 2/(Ra), 
their eigenfunctions can be chiral eigenstates, since D(A) and 75 commutes for these real 
modes. In the following, let us denote the number of the left(right)-handed zero eigenmodes 
as Nl(Nr) and that of the left(right)-handed A^- = 2/Ra (doubler) eigenmodes as n^nR). 

The propagator of the massive overlap fermion (for each flavor) can be, thus, expressed 
in terms of these eigenvalues and eigenfunctions as, 



Sa(x,v) = 

{n|ImA^>0} 



0n(^)0nd/) f , 750n(^)0n(?/) t 75 

f m X£ + m f m Xfi + m 



1\[A r ,A 

h rn + ti 2/(ite) ' 1 } 

where f m = l — Rma/2, and N R+L = N R + is the total number of zero modes while 
n R+L = n R + n L the total number of doubler modes. 

A measure for a given gauge field A can be also written in terms of eigenvalues as 

P m (A) = e- SY ^m N ^ + L(A R ) N f< + L J[ (Z 2 m X*X* + m 2 ) N f , (12) 

ImA„>0 

where Sym(A) is the gauge part of the action (whose explicit form is not needed in this 
work), A R = 2/(Ra) and = 1 — m 2 /A 2 R . Note that for even Nf, P m (A) is positive 
definite and an even function of m. 

It is important to note that all quantities which consists of Sa(x,v) and P m (A) are finite 
at V < 00, m 7^ and a / 0. We then carefully take the V — > 00 and m — > limits, 
without worrying about possible ultra-violet (UV) divergences, until we eventually take the 
continuum limit. 
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B. Chiral Ward-Takahashi identities on the lattice 

Now let us study the quantum aspects of the symmetry, performing the functional integral 
of an operator O over the quark fields, 

(0) F = J Oe~ SF . (13) 

The global lattice chiral rotation Eq. (J7J) gives the integrated Ward-Takahashi (WT) identity, 

((8 a0 J -8 a S F )O + 8 a O) F = 0, (14) 

where Jo is the contribution from the chiral anomaly, or the Jacobian of the measure, 

J = -2*iV 4 E E #v(s) f 75 (l " | od) <Pi(x) 

x N=n,k,K \ Z / 

= -2iN f x Q(A), (15) 

where Q(A) = N% — is the index of the overlap Dirac operator [14J], which gives an 
appropriate definition of the topological charge for the given gauge configuration A. 

In this paper, we consider the (volume-integrals of) scalar and pseudoscalar density op- 
erators 

S a = a A Y.W a {F{D{A))^}{x), (16) 

X 

P a = a 4 Yl$T a iK(F(D(AM(x), (17) 

X 

and their correlations. These two operators are transformed as 

5 b S a = 2 E d c ab P c , 5 b P a = -2 E d c ab S C) (18) 

c c 

where {T a ,T b } = 2^ c d c ab T c . In particular, in the Nf — 2 case, we have 

S b S a = 2S ab P , 5 b P a = -2(5 afe S , , (for a, b = 1, 2, 3), (19) 
8 S a = 5 a S = 2P a , 5 P a = 5 a P = -2S a , (for a = 0, 1, 2, 3), (20) 

where we adopted the normalization (T a ) 2 = 1 2X 2 without summation on a. It is now 
obvious that our mass term in the action Eq. (j3J) can be simply expressed by mSo, and its 
transformation is S a SF = 2mP a . 
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C. Basic properties and assumptions 



In this subsection, we explicitly give the basic properties and assumptions used in this 
paper. 

If the SU{2) L x SU{2) R chiral symmetry is restored at T > T c , we should have 

lim lim (S a O) m = (for a ^ 0), (21) 

m— >0 V— >oo 

for an arbitrary operator O, where an average over gauge fields is defined by 

(0(A)) m = 1 1 VAP m (A) O(A), Z = J VAP m (A). (22) 

Here we have put the subscript m to remind the readers of the m-dependence. 

In the following analysis, we will normalize the operator O (by multiplying 1/V with an 
integer k ) so that limv/^. 00 (5 a C) is well-defined. Note that P m (A) is positive for even Nf 
and JVAP rn (A)/Z = 1. 

In our analysis, we assume that the vacuum expectation values of m-independent observ- 
able O(A) is an analytic function of m 2 , if the chiral symmetry is restored. Therefore if 
O(A) is m-independent and positive for all A, and is shown to satisfy 

lim ±-(0(A) l °) m = (23) 

with a non- negative integer k and a positive integer Z 0; we can write 

(O(A) l0 ) m = m 2 ([ fe / 2 ] +1 ) J VAP(m 2 ,A)0(A) l \ (24) 

where [c] is the largest integer not larger than c, P(0, A) ^ for 3 A and 
J VAP(m 2 , A)0{A) l ° is non-negative and assumed to be finite in the large volume limit. In 
other words, the leading m dependence arises from the contribution of configurations which 
satisfy P(0, A) ± 0. 

Under the above assumption, it is easy to see 

(0(A) l ) m = m 2 ^' 2 ^ J VAP(m 2 ,A)0(A) 1 = 0(m 2 ™ +1 ^, (25) 

for an arbitrary positive integer /, as long as J VAP(m 2 , A)O(A) 1 is finite, since 0(A) l ° and 
O(A) 1 are both positive and therefore share the same support in the configuration space. 
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More generally, if a set of non-negative m-independent functions Gi(A) satisfies 
(Oi(A)) m = 0{m 2ni ) with non-negative integers rii (i — 1, 2, 3, • ■ • A;), it is easy to see that 

(f[O t (A)\ =0(m 2n — ), (26) 

where n max = max(ni, n 2 , • • • , n k ). 

If a non-negative Oo and an arbitrary operator 0\ are m-independent and satisfy 

(Oo) m = 0(m 2 -), (0 1 ) m = 0(m 2 -), (27) 

we then have 

(O O 1 ) m = m 2n ° J ' VAO Q {A)Oy{A) {P + (m 2 ,A) + P^(m 2 ,A)} = 0(m 2n °), (28) 

irrespective of values of no and rii, where P(m 2 , A) = P + (m 2 , A) + P_(m 2 , A) and 

1 (A)P + (m 2 ,A), 1 (A) > 
Oi(A)P_(m 2 ,A), Oi (A) < . (29) 
0, O 1 (A) = 



1 (A)P(m 2 ,A) 



As will be seen later, we have 



lim lim -^-{N^, T ) m = 0, (30) 

m->0V^oo mV R+L/m ' v ; 

as a constraint from the chiral symmetry restoration. This leads to 

lim hN* +L ) m = 0(m 2 ). (31) 

V— >oo V 

This condition is, however, much weaker than the naive expectation that the configuration 
A which give N^ +L = 0(V) has the weight P m (A) oc m N f°^ and therefore is much more 
suppressed in the large volume limit. We do not assume such a highly suppressed weight 
P m (A) in this paper. As will be shown later, however, we can further prove that 

lim ^« +L ) m = 0, (32) 

for small enough m, using our weaker assumption, Eq. (|24l) . 

Note that analyticity in m 2 for physical observables and its consequence Eq. (l2~4j) do not 
hold at T < T c , where the chiral symmetry is spontaneously broken. For example, the 
topological charge Q(A) is expected to satisfy 

lim hQ{A) 2 ) m = r ^ + 0{m% (33) 



where £ is the chiral condensate. The odd power of m reflects the non-analyticity of the 
QCD partition function at m — 0. 

In the following analysis, the thermodynamical limit of the eigenvalue density for a given 
configuration A, 

1 



p A (X) = lim - £ *(A-V^A£), (34) 

n(ImA£>0) 

plays a crucial role. Since the temperature of the system is fully controlled by P m (A), 
the eigenvalue density itself is not sensitive to the temperature. It is also notable that 
J^ R dX p(X) is finite on the lattice. Therefore, p A (X) is positive semi-definite for arbitrary 
choice of A and A. 1 

Although the original eigenvalue spectrum at finite V is sum of delta functions, we expect 
that such a spiky feature is smeared out in the thermodynamical limit, and p A (X) becomes 
a smooth function. We here further assume that p A (X) can be analytically expanded around 
A = 0, 2 

oo \ n 

P A W = E Pn^T- (35) 

An arbitrarily small convergence radius of this expansion, denoted by e, works well for 
our later discussion as we take the massless limit. As is well-known as the Banks-Casher 
relation [l5| and will be seen later, \im m ^ (p A } m 7^ means the spontaneous chiral symmetry 
breaking. 



III. CONSTRAINTS FROM THE SU(2) L x SU(2) R RESTORATION 

In the following analysis, we concentrate on the case with Nf = 2. In this section, we 
derive the constraints on the eigenvalue density of the Dirac operator in the SU (2) l x SU(2)r 
chiral symmetric phase at finite temperature. 

1 Strictly speaking, p A (X) has a logarithmic divergence in the continuum limit, which can be absorbed by 
multiplying the quark mass m. 

2 More precisely we here assume that configurations which do not have the expansion (|55j) are measure zero 
in the path integral with P m (A). This assumption excludes a possibility that (p A (X)} m has a fractional 
power such that (p A (X)) m ~ A 7 with non-integer 7 at small A. We consider the fractional case later in 
Sec. ED 
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A. WT identities for scalar and pseudo-scalar operators 



Let us consider a product of scalar and pseudoscalar operators defined in Eqs. (1161) and 

(n _ pn.1 C"i2 p"3 qri4 I Ofi\ 

^ni,n 2 ,n 3 ,n4 _ r a °a r Q °0 i \ ov J 

where a represents a non-singlet index (a = 1, 2, 3). Here and in the following, a summation 
over a is not taken, and we explicitly use "0" for the singlet operators. 
Non-trivial WT identities are obtained from a set 

q{n) = {O nun2in3j7l4 \ nx +n 2 = odd, m + n 3 = odd, = A}. (37) 

i 

which requires the operator to be non-singlet, and parity odd. More explicitly, we have at 
T>T C 

lim lim ^-(S a O nun2)n3! n 4 ) m = for 0„ 1 ,n 2 ,n 3 ,n4 e (38) 

where 

'in^(D ni ^ n2 j r \^ n3 —\^ nA + 2?T.4(!? ni ^i jn2jn3jn4 _x. (39) 

Note that the fermion integrals are performed before the gauge integrals: (C) m = ((0)p) m 
but we have omitted (■ ■ -)f for notational simplicity. Here the minimum power fc which 
make the V — >■ oo limit finite depends on the choice of O nijn2>n3>ni . See further details such 
as a relation of 7ii,7i2, 713,714 to in appendix I A 51 



B. Constraints at iV = 1 

At N — 1, there is only one operator Oiooo = -Pa in Ojf =1 \ which gives 

5 a P a = -2S . (40) 



Using the decomposition in Eq. (ITT]) , and normalization conditions that 
(75^, 750^) = ^nm, the tliermodynamical limit of the functional integral for So is expressed 
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as 



hm -<-S„>„ - hm - ^— + ^ ZIXK + ^ [I - -jj-j ^ 

where 

h = m £' d X A X)^^, *<«) = 1-^. (42) 

In the chiral limit that m — > 0, the only very vicinity of A = contributes to the integral, 
since 

W(A) (43) 

is finite for arbitrarily small but positive e, and thus, does not contribute to I\ in the limit. 
Expanding p A (X) for A < e, (see Eq. (155])). it is not difficult to obtain (see appendix I A 3p 

JO ^m* + m 

= irp A + 0(m). (44) 

As an exercise, let us consider the T < T c case, where the chiral symmetry is spon- 
taneously broken. Assuming that \\mv^oo{N^ +L ) m /V — >■ 3 , the famous Banks-Casher 



relation [15J is reproduced: 

lim Jim 1 r^r 7 (-S ) m = vr lim (p A ) m (= vr lim (p A {0)) m ) + 0. (45) 

On the other hand, in the chiral symmetric phase T > T c , we require 

lim lim ^(-S ) m = lim lim — ^(N£ , 7 ) m + JV/ hm = 0. (46) 

Since both of N^ +L and ii are positive, it is equivalent to separately require the two con- 
straints: 

hm ^f(N A +L ) m = 0(m 2 ), (p A } m = 0(m*). (47) 

Using Eqs. ( 126]) and (125]) (po)m = 0{m 2 ) implies (If) m = 0(m 2 ), which will be useful in 
the analysis below. 



3 In chiral perturbation theory, one can confirm (\Q(A)\) m /V is an 0(l/y/V) quantity (even when m is 
finite). Since the minimum of N^ +L is equal to |Q(A)| in the topological sector of Q(A), it is natural 
to assume (N^ +L ) m /V is also 0(1/W). Moreover, using the fact that there is no massless pole in the 
non-singlet scalar correlator (S a S a ) in the chiral limit, one can show that (Ng +L ) m /V = 0(m 2 ). 
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C. Contribution from zero modes at general N 



Before extending our analysis to higher N, let us discuss the fate of the zero-mode con- 
tribution at general N. For this purpose we consider an operator 0i,o,o,jv-i £ ^ a N \ whose 
non-singlet chiral WT identity requires 

lim lim (-(00,0,0,*)™ + ( N ~ 1) (02,o,o,tf-2)m) = 0. (48) 

in— >0 V— >>oo 



Its dominant contribution at large volume is 
and, therefore, from the positivity of N^ +L and 7i, 



(49) 



lim ((^A-L) Af )m = 0(m N + 2 ) (for even iV) 
y ™ ^ [0(m N+1 ) (for odd iV) 

Since this holds for arbitrary N, and N^ +L does not explicitly depend on m, we conclude 
that 

lim ( N R+ L )m = (51) 

at small but non-zero m. 

This result implies that any zero-mode's contributions to an arbitrary local operator are 
measure-zero in the thermodynamical limit, as we have already seen an example in Section U 
17| . Therefore, we hereafter set \uay^ f00 {N^ +L ) m /V = even at small but non-zero m. 



D. Constraints at N = 2 

We next consider the N = 2 case. In this case, two WT identities from 0iooi and 
Coiio £ 0^ =1 ) require that the so-called (non-singlet) chiral susceptibilities, 

X^ = j- 2 (S 2 - P 2 a ) m , X v ~ 8 = ^ ~ Sl) m (52) 

vanish in the V — > oo and m — > limits at T > T c . The first one, x"~ 11 ~> has been already 
examined in the previous subsection. 

In a similar way as the N = 1 case, x v ~ S can be expressed in terms of eigenvalues as 

lim X v ~ & = Jim (-^%Q(A) 2 ) + Nf ((^ + I 2 )\ , (53) 

13 



where h is defined by 



7, ^ ?' iX Ax) ^W^0 



7 + I)^ + ( 2 + ^- lo 4)^ + °W- < 54 > 



Noting 

h + I 2 = 2ml 3 = 4m 2 f dX /(A) f X * ] = + 2pf + 0(m) (55) 



m Jo (Z^X 2 + m 2 ) 2 

for an arbitrarily small (positive) parameter e, and expanding p A (X) around A = 0, we obtain 
a condition in the chiral limit that 



lim lim x v 5 = Nt lim 

m— ¥0 V^roo m— S-0 



lim Nf{QW U + £ <rf> + 2(pf)m + 0(m) 
y^oo 777/1/ m 



(56) 



Since we have already known (p A ) m — 0(m 2 ), this condition leads to 



lim N f (Q(Ap m = 2{p A )m + 0{m 2 ) (5?) 
y-s>oo m z V 

Note that there should be no 0(m) term in Eq. ( l57l) according to analyticity in m 2 . Therefore 
O(m) term can not be cancelled in Eq. ( 1561) at non-zero m. 



E. Constraints at N = 3 

From the WT identities at N = 3 except the one considered in the subsection IIII C\ the 
following four quantities, 

Xi = = «? </i/ 2 > m + 0(1/V), 

x2 = <o^ = _ iVf flW _ 2/ \ + 0(1/n 

\ I m 

\2\ \ 



(Co02l)r 



X4= (o^k = _^ (/f>m + 0(1/n (58) 

should vanish after taking the V — > 00 and m — )• limits. Here 7*3 (and its asymptotic form 
near the chiral limit) is given by 
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Substituting the explicit form of I± and I2, the result (1471) in the previous subsection and 
our assumptions in Eq. fl26|) and ff28l) give 

(hh) m = 0(m), {{hf) m = 0(m 2 ). (60) 

so that Xi an d Xi automatically vanish in the V — > 00 and m — > limits. Using the same 
assumptions and the result Eq. f[5"Tj) the following relations can also be shown 

«m). =I ww). +OMi (61) 

777/1/ m^l/ 

(QWm = 2 K>m + 0(m) (62) 



The remaining non-trivial two conditions are 

'(Po)m , (A 



lim lim X2 — —nNr lim 

m— >0 V— >oo m— >0 



+ 



m 2 2 



0, (63) 



lim lim X s = nNi lim lim v ^ v OT " = 0. (64) 

m->0 y->oo ' m->0 V->oo 171 V 

From the first condition, we obtain a constraint, 



(tf>m = -m 2 ^ + 0(m 4 ). (65) 

Moreover, since (pq ) m is positive (which is required by the positivity of (p A (0)) m ), (p^) m 
must be negative for small m. 
The condition for xz leads to 

lim {Q{A) y )m =0{m 2 ). (66) 

This condition does not necessarily give stronger constraint than {Q(A) 2 ) m = 0(m 2 V) and 
(Po )m — 0(m 2 ), since it only requires that a set of gauge configurations which satisfies both 
Q{Af ^ and p£ ^ has a weight m 4 P(A,m 2 ) + 0(m 6 ). 



F. Constraints at iV = 4 

The 8 WT identities at N = 4 give 7 independent constraints 

(04000 — C*0004)m — > 0, (04000 — 302OO2)m — ^ 0, 
(00400 — 0OO4o)m ~ ^ 0, (00400 ~~ 30o22o)7n ~~ ^ 0, 
(^2020 — 0O2O2)m ~ > 0, (02200 ~ 0OO22)m — > 0, 

(20 lm - 0Q2O2 + 0OO22)m 0, (67) 
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where the V — > oo and m — > limits are abbreviated by the arrows. The 0(y 4 ) contribution 
from Sq in the first equation has already been considered. 
At 0(V 3 ), there are 3 conditions, 

1 /c>2o2\ 7^3(^1)' 



lim ^(S 2 a S 2 ) m = -m{llh) m 0, (69) 



H-^oV^)^^-^}) -0. (70) 
v-kx> V/ 01 J \ [m m <rv J / m 

It is not difficult to confirm all of them are automatically satisfied, since 

(I?) m = ({vrp^ + 0(m)} n ) m = 0(m 2 ), (71) 
for any integer n > 2 from (po) m — 0(m 2 ), and 

(ggm» = ({^ + o(m)} 2 g(A) 2 ) m = 

from Eq. fl66|) together with the assumption (|28p for remaining cross terms. Namely, these 
three give no additional constraint. 
At 0(V 2 ) we have 

y^i^t ~ Pq)™ -> 0, — (S 4 - 3S' 2 P 2 ). m -> 0, — (P 2 {Pq - Si) - 2P a S a P S ) m -> 0. 

(73) 

After a little algebra using the formulas in appendices IA2I and IA 5[ the first condition 
becomes 

6jV 3 

3Nj((I 2 + h/m)(I 2 - h/m)) m + —L(Q(A) 2 h} m - -^L-(Q(Af) m 0. (74) 

Using 



6 " I = "» {-7n + 7 + ^) + + S ^ 4 ' 08 i 1 + 0(1) ' (?5) 



Eq. (I55p and ((Po) n ) m — 0(m 2 ), we can show that the first term in Eq. (174|) is at most loga- 
rithmically divergent in the limit m — > 0. Note that the second term is also logarithmically 
divergent due to cross contributions from the 0(m) terms in I\ (appendix IA 31) and Eq. (1571) . 
Therefore, in order to satisfy Eq. (1741) . the last term should not be power divergent and at 
least fulfill 

lim -^(Q(A) 4 > m = 0(m 4 ), (76) 
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which leads to 

^y-M A T k )™ = 0(m% (77) 

for an arbitrary positive integer k. Combining this with Eq. (1571) . we obtain a constraint on 
the spectral density, 

(pf) m = 0(m 2 ), (78) 

so that Eq. (1741) now becomes 

7T 2 TV 4 
- * N f-M) 2 )™ - 0. (79) 

Since both terms are negative semi-definite, this WT identity requires 

((po) k )m = 0(m 4 ), limi(Q(i) 2, ) m = 0(m 6 ), (80) 

V — 5>oo V 1 

for arbitrary positive integers k and /. The first condition also gives 

(p 2 A ) m = 0(m 2 ), (81) 

from Eq. (165|) . 

The last constraint Eq. (IHTj) can be obtained through a different argument. From Eq. (1781 . 
the eigenvalues density near chiral limit becomes 

(p A (A)) m = (p^y + 0(A 3 ) + 0(m 2 ). (82) 

The positivity of {p A (X)) m implies (p^) m > near m = but this contradicts with the 
positivity of (pg )m m Eq. ( 1651) . unless (p^) m = 0(m 2 ), and thus, (pg ) m = 0(m 4 ). 

It is now easy to see that the second and third conditions in Eq. ( 1731) are automatically 
satisfied: the second one gives 

3N 3 

6N 2 f m(I 2 I 3 ) m - -^-AhQ(Af) m = 0(m 2 ) + 0(m 4 ), (83) 

while the third one is evaluated as 

3iV? 

QN 2 f (hh) m - -^(hQ(A) 2 ) m = 0(m 2 ) + 0(m 3 ). (84) 
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G. Special constraints at general N 



In this subsection, we consider a special type of operators : 0o,i,(4fc-i),o £ & a N ^ & t a 
general positive integer k, whose non-singlet WT identity gives a condition, 



lim MAk - l)SiP^ - Pt) m = 0. (85) 

m— >0 



At the leading order of V (V 2k in this case), the above condition corresponds to 

- ^ ~ (hP£ k - 2 ) m - ^- k (Pt) m 0, (86) 



in the chiral limit. From the results in appendices IA 21 and IA 5\ 

7 4fc\ 2k /_ AT2n( A\2\ 2k ~ n 



V \ m?V I \ m 

where we used a definition (—1)!! = 1. The non-singlet WT identity is expressed by 

-,«- 1 ,E^«»-..: ! (m^(?)V fIl ) 



n=0 



n=0 



From the above condition(s), we would like to inductively prove that 

^ = 0(m 4fc + 2 ), ((p^> m = 0(m 2fc+2 ), (89) 

holds for an arbitrary positive integer I and k. 
Suppose 

^ = OK1, <(p^> m = 0(m 2fc ), (90) 

is obtained from the WTI at N = 4k — 4 (this is true for k = 2). The constraint above is 
then reduced to 



-A 2 g(A) 2N 

m 2 ^ 



(4fc-l) Hn^H N f {I 2 + 2kh/m)) ->0, (91) 



where only those terms with n = 0, 2/c — 1 in the first summation, n = 0, 1, 2/c in the second 
summation remain. While the first and the third term is finite and linearly divergent in the 
m — > limit, the second term is seen to be quadratically divergent from eq. (190]) as 

Nf (Q(A) 4 % 



m 4k y2k 
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O (m~ 2 ) . (92) 



In order that the WTI to be satisfied, the quadratic divergence should be absent, so that 
(Q(A) 21 /V l ) m = 0(m Ak ) for arbitrary positive integer I, thanks to Eq. ([25]) . 

Using this result the third term disappears faster than others and the WTI becomes 

_ (tt _ 1)!wri «^_5«^^ . (93) 

Note here that the both terms are negative semi-definite and therefore each term must vanish 
in the chiral limit. This completes the proof for Eq. (I59"j) . 

Since k can be arbitrarily large, we now have another stronger constraint on the zero- 
mode's contribution: 

lim = 0, (94) 

and that on the spectral density, 

(Po >m = 0, (95) 
which hold even at small but non-zero m. 



H. Short summary of the constraints 

Here we summarize the constraints obtained in this section. For the eigenvalue density, 
we have 

(f£) m = 0, (pf) m = 0(m 2 ), (p*) m = 0(m 2 ), (96) 

at a small but non-zero m. Namely, the eigenvalue density must have a form 

lim(/(A)) m = (pi) ^ + O(X% (97) 

We believe that this new condition is not only stronger than those found in the previous 
works, but also the strongest since we know that the Nf = 2 massless free quark theory has 
(p£)o keeping the exact chiral SU(2) L x SU(2) R (and U(1)a) symmetry. Therefore, 
it is very likely that we will not find any additional information from N > 5 correlation 
functions. 

For the discrete zero modes, we have obtained 

hm l J -((N£ +L ) k ) m = 0, hm -L(Q(Af k ) m = 0, (98) 

for an arbitrary positive integer k at a small but non-zero m. These zero-modes gives no 
contribution to the correlation functions we are considering. 
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IV. FATE OF U(1) A ANOMALY 



In this section, we discuss how the constraints in the SU(2)l x SU(2)r symmetric phase, 
obtained in the previous section, affect the U(1)a breaking correlators. Here we consider a 
set of (pseudoscalar singlet) operators, 

O ( N) = {O nun2jn3tni \m + n 2 = even, m + n 3 = odd, ^2m = N}, (99) 

i 

and its chiral U(1)a rotation, 

"0^ni,n2,M3,n4 '^'^'I^ni—l,n2+l,n3,n4 2?7<2^ni+l,n2 — 1,713,714 

2n3(9 nijn2jn3 _i n4 _(_i + 2?^40 r j, 1)r j, 2)n3 -|-i )n4 _i. (100) 
For later convenience, let us also define a set of (scalar singlet) operators, 

q{n) = {C ni)n2i „3 in4 1 ni + n 2 = even, ni + n 3 = even, ^ri; = iV}. (101) 

i 

As QCD keeps the vector like SU(2)v symmetry and the parity symmetry, any operator 
having non-zero expectation value must be a member of 0^ N >. Note that we have already 
introduced the set of pseudoscalar non-sing let operators { a N) in Eq. (157jl . 

Since the U(1)a transformation is anomalous, the expectation value of the variation 
(SoO) m is non-zero as the WT shows (see Eq. (TT5l)). 

\im (5 O) m = 2iN f ]im(Q(A)0) m , (102) 

771—5-0 771— >0 

and the U(1)a symmetry is broken. 

It is, however, still possible to have zeros on the both sides of Eq. f ll02p . If this is the 
case, the U(1)a anomaly is invisible. In fact, we show below that the constraints obtained 
in the previous section are strong enough to suppress the variation (S O) m for O G Oq N \ to 
be zero in the large volume V — > oo and chiral m — > limits. Namely, the U(1)a symmetry 
must be restored, at least, for the operator set Oq N \ 

A. Odd N case 

For the odd iV cases, we can show a relation for the number of operators \ \ = \Oq N ^ \ = 
\0^\, where \A\ denotes the number of independent operators in A. See appendix IB1 for 
the details. 
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As the exact chiral 577(2) l x SU (2) r symmetry requires | independent WT identities 
in the chiral limit to be zero, 

|0(iV)| 

lim(<M) m = lim MijiOjU = 0, { G Of >, 3 G (9 (7V) , (103) 

where the matrix M is specified by the WTI one considers. Since the chiral transformation 
keeps independency of the operator, it can be proved that det M ^ 4 , and the WTI requires 

lim (Oi) m = for any O t G {N) , (104) 

or equivalently, there is no operator in 0( N \5 a Of \ and 5 O { Q N) which has a non-zero 
expectation value in the SU(2) L x SU{2) R symmetric phase. 

Since the U(\)a variation of any operator in Oq N ^ is an element of 0^ N \ we can conclude 

hm(5 O) m = for any O G . (105) 

Without referring any specific constraints obtained in the previous section, we can thus show 
that the U(1)a breaking is invisible for these operators. 



B. TV = 2, 4, and 6 

For even N, the situation is not so simple as the odd TV's (see appendix [B] for the details.). 
We need to examine the WT identities explicitly. 

At N = 2, there remains one non-trivial susceptibility, but one can immediately show 
that it should vanish: 

1 Nf 
X*~ v = T/ {Pl ~ PoU = Jim -^-(Q{Af) m = 0, (106) 

for small but non-zero m, thanks to Eq. ((23). Noting P 2 - P 2 = (P 2 - Sg) + (Sg - S* 2 ) + 
(5 2 — P 2 ), we can also show 

X 5 - = ^- 2 (S 2 a - S 2 ) m = 0(1/V) + 0(m 2 ). (107) 

4 If dct M = 0, we can construct a chiral invariant operator from a linear combination of operators in 0^ N \ 
Since all operators in Oi N ^ have odd numbers of the index a, however, no chiral invariant operator should 
exist in Oa . 
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Since LHS of Eq. HlOTf) is the (double) volume average of LHS of Eq. ([T]), this give another 



proof that the U(1)a breaking effect in Ref. 17| cannot survive in the thermodynamical 
limit. 

At N = 4, there are two non-trivial susceptibilities, 

X5 = (^0022 — C2002)m, X& = (^0022 ~ Co22o)m- (108) 

Neglecting N^ +L /V and Q(A) 2 /V terms and using the constraint on I\ obtained in the 
previous section, both of them disappear as 

lim lim ^ = - lim lim N 3 ( ( + iX) = 0, (109) 
m-»o y->oo V 3 ™->o y->oo J \ m 2 V \ mV I 

Um Um " = u m u m ^ / (ULl + h ) 2 + h " Mm \ 



in 



->o y-s>oo V 3 m->o v-hx> m \\ mV j \ mV mV I / 

N 3 

= lim (if) = 0. (110) 
At iV = 6, we have 4 non-trivial susceptibilities, 

X7 = (Co024 _ C ) 2004)m, X8 = (^2040 — @0204:)m, (HI) 
X9 — (^0420 — 0OO42)m, XlO = (00042 _ Co024)m- (H2) 

In the large volume limit, they behave as 



lim — 

V->oo V^ 5 


= o, 




(113) 


Iim — 


= N 5 f (I 2 ll) m 


= 0(m 4 ), 


(114) 


Iim — 

y^oo ^/4 


Ni 


= 0(m 2 ), 


(115) 


r Xio 
lim — - 


N f 5 

m 


= 0(m 4 ), 


(116) 



all of which vanish after the chiral limit is taken. 

We thus conclude that the U(1)a symmetry breaking is not viable at least for N < 6. 

C. General even N 

In order to consider the general iV case, let us look at RHS of Eq. f ll02p . Namely if we 
can show 

limlim ±(Q{A)O) m = 0, (117) 

m— >0 V— >oo v 
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with some appropriate power of k, we can prove that LHS of Eq. fll02j) also vanishes. In the 
analysis below, we divide Oq into two classes : the one with (m, 712, 7i3)=(even, even, odd), 
and another with (rii, ri 2 , ri3) = (odd, odd, even). 

For the former class, or more explicitly in the case of (711,712,713) = (2/c 1; 2/c 2 , 2/c 3 + 
1), the leading contribution in S O nitn2>n3>ni comes from — n 3 O nitn2>n3 ^i jTlA+ i, whose leading 
contribution in V has 0{V k ) with k = k\ + & 2 + k% + + 1 (see appendix IA 5 j) . Therefore, 
we have 

l ^±(Q{A)O niMni ) m ~ iV / fc+1 r i3 (2A; 1 -l)!!(2A ;2 -l)!!(2A;3-l)!! 

Q2 sj^h+ks 



*V^^J <-'»)"■ <-'>>"7 ins) 

in the large volume limit, where zero modes contributions are neglected except the first 
term. According to the property Eq. ( 1981) and assumption ( 1281) . the RHS indeed vanishes 
inl/^oo limit at small but non-zero m. In the case with (m, n 2 , n 3 ) = (odd, odd, even), a 
similar analysis gives the same conclusion. 

We conclude that, for a class of operators we have considered in this paper, the U(1)a 
breaking effects are invisible in the thermodynamical limit. 



V. POSSIBLE ARTIFACTS IN LATTICE QCD 

In the previous sections, we have investigated the symmetry restoration for T > T c , fully 
relying on the exact chiral SU(2)l x SU(2)r symmetry and taking the thermodynamical 
limit V — > 00. However, numerical lattice QCD simulations must be performed on a finite 
volume, sometimes employing fermion actions which explicitly break the chiral symmetry. 
In this section, we would like to briefly address possible systematic effects of not having 
these two key properties. 

First, we discuss the explicit breaking of the chiral SU{2)l x SU(2)r symmetry. In order 
to characterize its violation, let us introduce a mass parameter mb rca k- Since mb re ak should 
disappear in the continuum limit, it is natural to assume mbreak ~ Aqcd ' where Aqcd 
is the QCD scale. At finite temperature T, we also have a possibility of mbreak ~ T 2 a or 
TTibreak ~ Aqcd^o- But it is unlikely that the lattice artifacts grow with T, since they are 
naively expected to be milder in the weakly coupled region at higher temperature. We also 
neglect the possibility of m bre ak ~ T% a since T c is not essentially different from Aq C d- 
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If one employs improved Wilson-type actions or staggered-type actions, it may be reduced 
to 0(a 2 ): mbreak ~ A.Q CD a 2 . F° r the domain- wall fermion action, as an approximation of 
the overlap fermion action, the suppression of the discretization effects could be stronger. In 
this case, the so-called residual mass, m Tes , is a good estimate for mb rea k- For this reason, we 
have introduced a rather abstract parameter mb rea k; to treat the conditions with different 
actions in a uniform manner. 

Now the discussion is simple. As loosing the required exact chiral symmetry, every result 
in the previous sections should be, in principle, modified by the effects of mbreak, unless 
some special cancellation mechanism occurs. The condition, {po) m = 0, could be an only 
exception, being as a definition of the symmetry restoration. 

Namely, instead of Eqs. (J57|) . we should have 

A 3 

(p A (X)) m = am brC akAQCDA + /3m br0 akA 2 + ((P3 )o + 7 m break/A QC D) ^j" H , (119) 

where a,j3, 7, ••• are unknown dimensionless 0(1) coefficients. Similarly, (Njt + i,) m /V, 
(Q(A) 2 ) m /V , x n ~ v an d so on should not be zero but 0(mb rea k)- Note here there is no 
reason for the U(1)a symmetry to get restored as long as the chiral symmetry is explicitly 
broken by the quark action at finite lattice spacings. 

Next, let us discuss the finite volume effects. Even in the continuum limit, it is possible 
to obtain different results from what we have shown in this paper. It is good to remember 
that Eq. (JTJ is an example of the U(1)a breaking in the continuum theory, as a finite volume 
effect. The recovery of the U(1)a symmetry is not as strong as the other symmetries, under 
which the Lagrangian is strictly constrained, but rather a consequence from QCD dynamics, 
which is manifest only in the thermodynamical limit. 

Above the critical temperature, the long-range physics of the system would be character- 
ized by the correlation length £, or inverse of some screening mass, which diverges at T c if 
the phase transition is of second order while remains finite for the first order transition. 

It is then natural to assume that a lattice QCD simulation has finite volume effects as 
functions of ^/V 1 / 4 . It is also important to note that, unlike the truly local quantity, whose 
finite volume effects are exponential ~ exp(— V 1 ^ 4 /^), the susceptibilities considered in this 
paper are volume averaged quantities, so that any finite volume effect is expected to be a 
power function of ^/V l / A . A careful estimate for the thermodynamical limit is thus required, 
in particular, for temperature near T c , where £ could become larger. 
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VI. EIGENVALUE DENSITY WITH FRACTIONAL POWER 



So far, we have assumed that p A (X) is analytic around A = 0, and used the expansion in 
Eq. f[3"5"j) . In this section, let us extend our analysis to a non-analytic case that 



p A (X) = c A X\ 



(120) 



with a fractional power 7 for A < e, where c A is a A-dependent constant. Since 
lim m _ i , (p A (0)) m = in the SU{2) L x SU{2) R symmetric phase, 7 should be positive as 
long as (c A ) m = 0(1). It is still true in this case that the only very vicinity of A = 
contributes to the WTI's. We thus can neglect additional terms with higher order fractional 
powers, even if they exist in the bulk region of A > e. 
In this case, Ii, I 2 , I3 are expressed as 



h 



— + I 2 = 2ml 3 = c A 
m 



m 7 (di + 0(m 2 )) + m (e x + 0(m 2 )) 

m?- 1 (d 2 + 0(m 2 )) + m 2 (e 2 + 0(m 2 )]) 



where d^s and e^'s are given by 



di = 7rsec 



77T 



' 77T \ 

d 2 = (1 -7)7rsec ( — J 



and 



ei 



e 2 



dXp A (X) 



2^o(A 5 



Z 2 A 2 + m 2 



+ e 



A ^Ap^(A) ^ 



r(i±i) 
3 (7 + 3)r to 



X 



(Z^A 2 + m 2 ) 2 
2 



+ e 



2( 7 -3)r(^, 

with the UV cut-off and the IR cut-off e. Note that d{S and e^'s are all finite. 



( 7 + 1)( 7 - 1) - 2-^(7 + 1)( 7 - 3) + ^(7 - 1)(7 - 3) 



;i2i) 
122) 



(123) 



(124) 



(125) 



A. < 7 < 1 



We first consider the case with 7 < 1. With the above expressions for /j's, let us reexamine 
the WTI's given in the previous sections. 
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For (Sq }m/y with an arbitrary N, the WT identity requires 



lim lim + \ =0, (126) 

m-X) V^oo \ mV 1 



III 



where the positivity of each term implies 



at small m, and 



lim {N£ +L ) m = 

V-s-oo V 



lim m N "<d?((c A ) N ) m = 0, (128) 

m— >0 



which is automatically satisfied for positive 7. 
At N = 2, we have 



■h +h .mf \ = ^ v>m _A&M=_> . (129) 

m m z v / m z v 



in 

Taking into account a fact that both c A and Q 2 are mass independent and their expectation 
value should be written as an even power of mass, both terms in Eq. HI 29 j) should vanish 
separately, which leads to 

(c A ) m = 0{m% (g( ^ )2)m = 0(m 4 ). (130) 

At N = 3, it is not difficult to see that all the non-trivial conditions are automatically 
satisfied with the above constraints. 

At N = 4 there remains one non-trivial WT identity: 

{6mW2 _ h/m))m + - «X ^ o. (131) 

Since the first two terms vanish in the chiral limit, we obtain a new constraint that 

lim (Q(A) 2 ) m = 6 2 

V— >oo v 

Let us finally consider the WT identity from Q^ =Ak = C ) o,i,(4fc-i),o as before. In a very 
similar way as the one in Section III I Gl we can show 

(c A ) m = 0, (134) 

even at non-zero m. This means that the Dirac eigenvalue density with a fractional power, 
Eq. (I120p . is incompatible with the SU(2)l x SU(2)r chiral symmetry restoration for < 
7 < 1. 
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B. 1< 7 < 2 



Next, let us consider 1 < 7 < 2. Our strategy is the same as in the previous subsection, 
except that the leading term is not 0(m 7 ) but 0(m) in 1^. 

Up to N = 4, one can easily confirm that most of the conditions are automatically 
satisfied for 1 < 7 < 2, keeping the constraints on the zero- mode contribution Eqs. (j!27p 
and (11331) unchanged. The only non-trivial WT identity appears at N = 3: 

lim^ = 2N f (I 3 ) m = 2N f d 3 m^ 2 (c A ) m + O(m) = 0, (135) 

which leads to a constraint 

(c A ) m = 0(m 2 ). (136) 
Namely, the fractional power 7 < 2 cannot survive in the chiral limit. 

C. 2 < 7 < 3 

In this case, all the conditions from WTI's are automatically satisfied up to N = 6 as 
long as Eqs. ( 1127ft and (11331) are satisfied. We thus have no constraint on (c A ). 

It is, however, important to note that excluding 7 < 2 in the chiral limit is enough to 
achieve all the U(1)a symmetric identities in Section HVl As discussed in subsection II V CI the 
zero-mode's contribution plays a more important role than bulk contributions from non-zero 
modes. 

VII. SUMMARY AND DISCUSSION 

In this paper, we have investigated the eigenvalue density p A (X) of the Dirac operator 
in the chiral SU(2)l x SU(2)r symmetric phase at finite temperature. In order to avoid 
possible ultra-violet divergences, we have worked analytically on a lattice, employing the 
overlap Dirac operator, which ensures the exact chiral symmetry at finite lattice spacings. 

From the various WT identities of the scalar and pseudoscalar operators, we have shown 
that a behavior such as (p A (\)) m oc A 7 for small A cannot survive in the chiral limit for 
7 < 2. If (p A (\)) m is analytical around A = 0, this means that it should start with a 
cubic term, as is the case with the free quark theory. Moreover, we have found a strong 
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suppression on the zero-mode's contributions in the thermodynamical limit. As shown in 
Eq. (]98p . they disappear even with small but finite m. It is worth mentioning that the 
use of the overlap fermion is crucial to obtain results in this paper since only this fermion 
formulation can preserve the exact chiral symmetry with non-perturbative cut-off, which 
makes our arguments more rigorous. 

The obtained constraints on the Dirac spectrum are strong enough for all the U(1)a 
breaking effects among correlation functions of scalar and pseudo scalar operators considered 
in this work to vanish in the limits of V — > oo and then m — y 0. Namely, there is no remnant 
of the U(1)a anomaly above the critical temperature at least in these correlation functions. 

This does not contradict with the apparently opposite results about the ^(l)^ restoration 
in the previous works. As we have shown in Section (H their U(1)a breaking parts cannot 
survive in the thermodynamical limit (V —> oo), but they could be finite on a finite box, 
which may be a part of difficulties for numerical lattice QCD simulations. 

We only use a part of chiral Ward-Takahshi identities to derive the constraints in this 
paper, which are therefore necessary conditions to be fulfilled if the chiral symmetry is 
restored. Our results strongly rely on analyticity in m 2 for m-independent observables 
and its consequence eq. (124)) . If our results are shown to be incorrect by some numerical 
simulations, these assumptions must be also violated in the simulations. 

One of the most important consequence of our study is that, since the U(\)a anomaly 
effect disappears in scalar and pseudo scalar sectors at T c , the chiral phase transition for 2 



flavor QCD is likely 



universality class 19 



;o be of first orderjlj or of second order in the U(2)l ® U(2)r/U(2) 



2 If . contrary to the expectation that the chiral phase transition of 2 



flavor QCD belongs to the 0(4) universality class. 
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Appendix A: Useful formulas 

1. Quark contractions for the (pseudo) scalar operator 

We give a contraction formula for the (pseudo) scalar operators when we integrate the 
fermion fields out. 

(S )f = -N f tvS A , (P ) F = -iN f tT"f 5 S A , 

(S 2 a ) F = -N f tr S% (P 2 ) F = N f tr ( l5 S A f , 

(S 2 ) F = -N f trS 2 A +(N f trS A ) 2 , (P 2 ) F = N f tr ( l5 S A ) 2 - (N f tr l5 S A )\ 

(S P ) F = -iN f tr l5 S 2 A + iN 2 f tvS A tv l5 S A , (S a P a ) F = -iN f tr l5 S 2 A , (Al) 

where S A (x,y) = F(D)S A (x,y). 

2. Trace of fermion propagators 

We give useful formulas for the trace of fermion propagators in a form of the eigenvalue 
decomposition. 

Let us first define S A as 



S A = J d A xid i X2 ■ ■ ■ d A x n S A (x 1: x 2 )S A {x 2 , x 3 ) ■■ ■ S A (x n , x x ) 

/n 
J[d 4 XiS A (xi,x i+1 ), (x n+1 =x 1 ). (A2) 



8=1 



Inserting the eigenvalue decomposition for the fermion propagator Eq. ffTTj) . we obtain 

1 ~ N A 1 ~ N A 

~trS A = + itrSi = % + / 2 , (A3) 

V mv V m z V 

I N A T 1 N A 

T^tr ( 75 ^) 2 = ^ + \ U^S A1 ,~S 2 A = ~^ + h, (A4) 

V m z V m V m 6 V 

-;tY-f 5 S A = — =, -j=tT^S A = (A5) 
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where U (i = 
limit as 



1,2,3) are expressed in terms of the eigenvalue density in the large volume 



hk-l 

So(A 2 ) 
h 







1 - 



A 2 ' 

ly R 



1 - 



m 



m 



-— + 2mJ 3 . 
m 



(A6) 

(A7) 
(A8) 



3. Various integrals of eigenvalue density 

The above Jj's are evaluated by expanding the eigenvalue density as Eq. f)35p . In evalu- 
ating -?2fc-i it may be better to rewrite hk-i = I^k-i + 0(m) where 

2^) 



T k-i - -^izi [ rf Ap A (A)— 

z m Jo (A 2 + m R ) 



by using = m/Z m . The expansion is given by 



(A9) 



re - V n A 7" (n) 

n=0 



with the expansion coefficient 



r(n) _J_ 



\ n 2m R g k {\^ 



dX 



n\ (A 2 + m|) fc ' 
At n > 2k — 1 the leading term in m is given by 

,n-2fc / ''^ 



T(n)leading 
'2fe-l 



2m 



R 



A' 



o n! 



A 2 



The explicit form of the coefficient is given as follows for a few k and n 

.2 \ 



(0) 



2 



r(l) m /? 

m R 



n 



(2) 



I" 



Afj tal1 



/ e \ _ m R e 
\mj K 2 R 



a! log U + 1 )-T? 



l R / 

3m| - e s 



3A| 



R 

( m %\ _i 

m fi II + -^2- I tan 



r(0) 



1 



Zm m R 



3 i) t 1 + S) tan " (i 



(-) 

\m R / 



m R e 
m R + e 2 



(A10) 



(All) 



(A12) 



(A13) 



(A14) 



(A15) 



2m R e 



(A16) 
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r(l) 



r(2) 



r(3) 




According to the equality the coefficient for I2 is given by 
2r 



r(0) 



r(l) 



r(2) 



(A|+m|) 2 



2m| 



A|) -6m i? (A| + m2 ? )tan- 1 (^) 



2(A 2 R +<) 2 
e 2 +m'j i 



+ 2m% + A%)- (A| + 4m£) (A| + m R ) logK + l 



3Af^ 



-e 3(A 



3m R (A^ + m^)(2A^ + 5m#) tan 1 ( 




4m R + m 



A1 



m 



«l2 



+ m 



-e 2 (A 



2m 



R) 



R 



(A17) 



(A18) 



(A19) 



(A20) 
, (A21) 



(A22) 



4. Integrals of eigenvalue density with fractional power 



If we consider the fractional power for the eigenvalue density with 7 > the eigenvalue 
integral 

2^(A 2 ) 



r(7) 
l 2k-l 



m R 



72fe-l 



dA A 7 - 



2 \fc 



'0 (A 2 + m|) 
is given in terms of the hyper geometric function as 

2 ' 2 ' m 



(A23) 



r(7) 
! 1 



r(7) 



2e 7 +M( 7 + 3) 2 F 1 (l,^±l-^- 



A| (7+1)2Fl ( 1 '^- ; ^- ; -m| / 



(7 + 1)(7 + 3)Z m m R 



(A24) 



e 7+l 
m R^m 



m. 



e 2 + m 2 



2\2 



7-1 



R 



2e 2 7 + 1 
A|7 + 3 



^1 1, 



7 + 1 
7 + 3 7 + 5 



iFi 1 



7+1 7+3 



ml 



e 4 7 



??2 



A* T + 5 



A 1 



7 + 5 7 + 7 



m 



(A25) 
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where 2 -Pi is the Gaussian hyper-geometric function given by 

2 Fi(a,/?,7;z) 



r( 7 ) ~ r(a + n)r(/3 + n)^ 



r(a)r(/9)^ r( 7 + n) n! 
Performing an expansion for m#/e < 1 we have 



I r ( 



n 



(7) _ . m R 



6 2 r(gr 



^r(2±i 



+ 



(=r) ! 



2^r(i- i )r( i + |) //m R y- 

(1 + 7)2™ IA e J , 



r(7) 



( 7 + 3)^r(2=i) 



4(7 - 3)^3 r (2±*) 



2e^ 



(A26) 



(A27) 



( 7 + 1)( 7 - 1) - -5.(7 + 1)( 7 - 3) + -r( 7 - 1)(7 - 3) 



+ o 



7-2 



V e / 



V e / 



:i-7)^ a r(V)r(ip) 



2Zi 



\ e / 



(A28) 



5. General correlation functions in large volume limit 

Here we consider the leading volume scaling of the general correlation functions made of 
SaS and P a 's. 

At given order N = 2(ki + k 2 + k^) + 714, there are two types of parity-even and chiral 
symmetric operators, 



Of = e> 2fc i+l,2/c2+l,2/c3+l,rt4— 3- 



For 0\ , the integration over the fermion fields in the large volume limit is given by 



(A29) 



"4 

F 



(A30) 



where each (O) p gives an 0(V) contribution. An overall constant coming from combinatorial 
factors is omitted here and hereafter. Therefore the leading contribution is 0(V kl+k2+k ' i+n4 ) 
for Of. 

Similarly we have 



(0?)f ~ (PaS^FiP^^iSD^iP^FiSoY, 



114—3 

F > 



(A31) 
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where 



(Po k3+l )F ~ (Po)F(P 2k:i )F - (Po)f(P^ f s . (A32) 

Since (P a S a )F and (Pq)f are proportional to Q(A) and is therefore 0(W), the leading 
volume dependence of Of is 0(V kl+k2+hi+n4 - 2 ). 

Appendix B: Structure of Ward-Takahashi identities for scalar and pseudo-scalar 
operators 

In this appendix, we summarize the general structures of WT identities among the scalar 
and pseudo-scalar operators, O nimm>ni = P™ 1 $a 2 Po 3 ^o 4 ■ 

In this work, we study the relation between three operator sets below, 

q{n) = {C niin , 2in3in4 | m + n 2 = odd, m + n 3 = odd, Y, n i = JV}, (Bl) 

i 

O { N) = {0 ni ,n 2 ,n 3 ,n 4 | ™i + n 2 = even, m + n 3 = odd, ^ m = N}, (B2) 
qW = {O nun ^ n:jyn4 \ni + n 2 = even, n x + n 3 = even, ^2m = N}. (B3) 

i 

Note that the only can have a non-zero expectation value in QCD with 2 degenerate 
quarks. 
Since 

i„OW = G O^} G (7V) , ^oO^ = {5 O|O G 0^} G (9 (7V) , (B4) 

our goal of this paper is to understand constraints from the SU(2)l x SU(2)r symmetry 
restoration, (<5 a (9f ) = 0, and to examine whether S O { N) can have a non-zero expectation 
value, which means that the U(1)a is still broken. For simplicity, hereafter we denote 
n = {nin 2 n^n^ instead of O ni>n2>n3>n4 to represent an operator. 

1. WTI's at odd N 

As shown in Section IIVI we can show {6 O ( N) } = if (6 a Oi N } ) = when N is odd. 
This follows from the fact that \0^\ = \O { N) \ = \0^\, where \0\ means a number of 
independent operators in O. Here we give a proof for this equality for general odd N = 2k+\. 
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At k = 0, we have only one operator for each set: = (1000) for 0^ N \ Ub — (0010) for 
O { N \ and n c = (0001) for O^. Thus, \0^\ = \O { N) \ = \0^\ = 1 for k = 0. 

At k = 1, we can create the operators by adding to the above %(X = A, B,C) a 
pair of the same operators, namely, adding 2 to one element of nx- Since each rix has 4 
elements, we have 4 operators for each set. We should, however, note that there is one 
additional type of operators for each set: ha = (0111) G 0[ N \ ub = (H01) G Oq N \ and 
nc = (1110) G C (Ar) . Therefore, \0^\ = \O { N) \ = \OW\ = 5 in total. For example, we 
have (1002), (1020), (1200), (3000), and (0111) in 0[ N \ 

In fact, every operator at higher k can be generated by adding 2 to one element of nx k 
times or adding 2 to one element of nx k — 1 times. A number of independent operators at 
a given k, therefore, is obtained by a number of selecting k (or k — 1) numbers from 1, 2, 3, 4, 
which is ^+3^3 (ah^Cs)- In total we have 

r , r (k + 2)(k + i)(2k + 3) 

fc+3^3 +k+2 t>3 — ^ 1 \Pb) 

for each set. This completes the proof for = | Cq | = at an arbitrary odd 

number N. 

2. WTI's at even N 

We next consider the case with N = 2k. We have nA{ns) = 0110(0011) and n^^s) = 
1001(1100) at k = 1 for OW (Of ] ), while n c = 0000 at k = and n c = 1111 at k = 2 for 
C?^) . As before, it is easy to count a number of independent operators for each case. There 
are 2 x (^+2^3) operators for 0f> and O ( N) , while k+3 C 3 + k+1 C 3 for O^l From this, it is 
easy to see 

\0 (N) \ - \Oj$\ = k+3 C 3 + k+1 C 3 - 2 x ( k+2 C 3 ) = (k+ 1) > 0, (B6) 

which means that > \0^\ = \O ( N) \ at N = 2k. Therefore 5 a O ( N) = is not 

equivalent to 6 Oq N ^ = 0. 
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3. Explicit WTI's at small N = 2k 



a. k = 1 

In this case two non-singlet WTI is given by 

(0020) - (0200) = 0, (2000) - (0002) = 0, (B7) 

while the singlet ones give 

5°(0011) = (0020) - (0002), (B8) 
5°(1100) = (2000) - (0200) = (0002) - (0020) = -5°(0011). (B9) 

Therefore, one non-trivial U(1)a rotation can remain. For simplicity, we omit the bracket 
of ({ni^n^n^)) here. 



In this case, 8 non-singlet WTI's read 



(2020) - (2200) - 2(1111) = 0, (4000) - 3(2002) = 0, (B10) 

(2200) - (0202) + 2(1111) = 0, 3(0220) - (0400) = 0, (Bll) 

(2020) - (0022) - 2(1111) = 0, (0040) - 3(0220) = 0, (B12) 

(0022) - (0202) + 2(1111) = 0, 3(2002) - (0004) = 0, (B13) 

which can be reduced to 

(4000) = (0004) = 3(2002), (0400) = (0040) = 3(0220), (B14) 
(2020) = (0202), (2200) = (0022), 2(1111) = (0202) - (0022). (B15) 

Note that these are symmetric under ri\ -H- 714, n 2 -H- n 3 . 

Using the above conditions, we have two independent quantities 

(0022) - (2002), (0022) - (0220), (B16) 

to examine the U(1)a chiral symmetry. 
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c. k = 3 



We have 20 WTI's from the non-singlet chiral symmetry: 

(4020) - (4200) - 4(3111) = 0, (6000) - 5(4002) = 0, (B17) 

(2400) - (0402) + 4(1311) = 0, (0600) - 5(0420) = 0, (B18) 

3(2220) - (2400) - 2(1311) = 0, (4200) - 3(2202) + 2(3111) = 0, (B19) 

(2040) - 3(2220) - 2(1131) = 0, (4020) - 3(2022) - 2(3111) = (B20) 

plus equations derived from the above by fiiOii4,n 2 Oii3, and 

(2022) - (2202) - 2(1113) + 2(3111) = 0, 3(4002) - 3(2004) = 0, (B21) 
(2220) -(0222) + 2(1131)- 2(1311) = 0, 3(0240) - 3(0420) = 0. (B22) 

The above conditions are summarized as 

(4002) = (2004), (0240) = (0420), (4020) = (0204), 

(4200) = (0024), (2400) = (0042), (0402) = (2040), 

(6000) = (0006) = 5(2004), (0600) = (0060) = 5(0420), 
6(2220) = 6(0222) = (0042) + (2040), 6(2202) = 6(2022) = (0204) + (0024), 
4(3111) = 4(1113) = (0204) - (0024), 4(1311) = 4(1131) = (2040) - (0042). 

(B23) 

In this case, there remain 4 non-trivial chiral U(1)a rotations as 

(0024) - (2004), (2040) - (0204), (0420) - (0042), (0042) - (0024). (B24) 
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